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ABSTRACT 

We calculate the probability to form giant radio halos (~1 Mpc size) as a function of the mass of the host 
clusters by using a Statistical Magneto- Turbulent Model (Cassano & Brunetti, these proceedings). We show that 
the expectations of this model are in good agreement with the observations for viable values of the parameters. In 
particular, the abrupt increase of the probability to find radio halos in the more massive galaxy clusters (M^2 x 
IO^^Mq) can be well reproduced. We calculate the evolution with redshift of such a probability and find that giant 
radio halos can be powered by particle acceleration due to MHD turbulence up to z~0.5 in a KCDM cosmology. 
Finally, we calculate the expected Luminosity Functions of radio halos (RHLFs). At variance with previous studies, 
the shape of our RHLFs is characterized by the presence of a cut-off at low synchrotron powers which reflects the 
inefficiency of particle acceleration in the case of less massive galaxy clusters. 
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I. Introduction 

Radio observations of galaxy clusters indicate that 
the detection rate of radio halos (RHs) shows an abrupt 
increase with increasing the X-ray luminosity of the 
host clusters. In particular about 30-35% of the galaxy 
clusters with X-ray luminosity larger than 10''^ erg/s 
show diffuse non-thermal radio emission (Giovannini 
& Feretti 2002); these clusters have also high temper- 
ature (kT > 7 keV) and large mass 2x IO^^Mq). 
Furthermore giant RHs are frequently found in merging 
clusters (e.g., Schuecker et al 2001). These observations 
suggest that there is a connection between thermal and 
non-thermal phenomena in galaxy clusters. 
Recent papers (Ensslin and Rottgering 2003; Kuo et 
al. 2004) have investigated the statistics of RHs and 
their connection with the thermal properties of the host 
clusters from a theoretical point of view. These works 
are based on assumpions in defining the condition of 
RHs formation from observational correlations and/or 
mass thresholds. Present data suggest that giant RHs 
may be accounted for by synchrotron emission from 
relativistic electrons reaccelerated by the turbulence 
generated in the cluster volume during merger events 
(Brunetti 2003; Brunetti this proceedings). Thus, with 
the aim to investigate the statistical properties and the 
connection between thermal and non-thermal phenom- 
ena in galaxy clusters, we have developed a statistical 
magneto-turbulent model (Cassano & Brunetti 2004, 
C&B model; Cassano & Brunetti these proceedings) 
in which we follow the formation of clusters of galaxies 
(making use of the extended Press & Schechter (1974) 
formalism) and estimate the injection of fluid turbu- 
lence and of fast magnetosonic (MS) waves during clus- 
ter mergers. Then we calculate the evolution of the 
electron spectra in the ICM and the resulting radio 



(synchrotron) and hard X-ray (Inverse Compton) emis- 
sion spectra. By using our model we can thus investi- 
gate the probability of formation of RHs in a well de- 
fined physical framework, the evolution with redshift of 
such a probability and the expected Luminosity Func- 
tions of RHs (RHLFs) . Here we applay the C&B model 
under the following assumptions: 

- We focus the attention on the formation of giant 
RHs only (radius Rh ^ 500 h^^ Kpc). 

- The magnetic field strength averaged over the 
emitting volume is assumed to be < _B >~ 0.5 
/iG, independent on the mass of the parent clus- 
ter. 

The adopted cosmologies are: EdS [Ho = 50 Km 
s-^Mpc-\ no,rn = 1) and ACDM {Ho = 70 Km 
s-Hdpc~\ no,„, = 0.3, = 0.7, as = 0.9). 

II. Occurence of RHs: predictions vs observa- 
tions 

By making use of the C&B model we have run Monte 
Carlo simulations to obtain a sufficiently large number 
of merger trees in order to have a large synthetic pop- 
ulation of galaxy clusters with a wide range of present 
day masses and temperatures. In this way we are 
able to statistically follow the cosmological evolution 
of the non-thermal emission and of the properties of 
the thermal ICM. Clusters with RHs in our synthetic 
population are identified with those objects with a syn- 
chrotron cut-off i/f,^10^ MHz in a region of 1 Mpc h^Q 
size. We have calculated the probability to form RHs 
(z < 0.2) in two mass bins: binA=[1.8 - 3.6] x IO^^Mq 
andbinB=[O.9-1.8]xlOi'^Af0 (EdS cosmology). These 
mass bins are consistent with those considerated in the 
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1. — Expected formation probability of RHs 
500ft.^Q^kpc, B 0.5/iG) as a function of param- 



Fig. 

{Rh 

eter rjt in a EdS cosmology (solid lines with error bars) 
and in a ACDM cosmology (dotted lines) in the mass 
bins: binA=[1.8 - 3.6] lO^^MQh^^ and binB=[0.9 - 
1.8] IOI^Mq/i^o^ for EdS case and binA=[1.9-3.8] •IQi'^ 
M0hfo^ and binB=[0.945 - 1.9] • lO^^ M^hj^ for the 
ACDM model. The two bottom dashed lines mark the 
observed probabilities for RHs in the mass binB while 
the two top dashed lines mark the observed probabil- 
ities in the mass binA; observational regions already 
account for Icr errors. 

observational studies and thus allow us to compare our 
expectations with observations. In Fig.(l) we report 
the probability to form a giant (~1 Mpc h^fl size) RH 
(red points) in the two mass-bins (including the statisti- 
cal error estimated from our Montecarlo simulations) as 
a function of the parameter rjt, which gives the fraction 
of energy of the turbulent motions injected by cluster 
merger which is channeled in the form of MS waves in 
the C&B model. 

The dashed blue lines mark the range of the observed 
probabilities (Giovannini et al. 1999) in the binA (top 
dashed region) and in the binB (bottom dashed region) , 
respectively. For a comparison in Fig.(l) we also report 
the probability to form a RH (z < 0.2) in a ACDM 
cosmology (green dotted hues). As expected, we find 
that at z < 0.2 the results are relatively indipendcnt 
from the considered cosmology, with the ACDM model 
being slightly less efficient. The main result is that 
in both EdS and ACDM models it is possible to find 
a unique interval of rjt in which the model reproduces 
the observed probability for both the cluster-mass bins. 
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Fig. 2. — Probability to form giant RHs (> IMpc h^Q 
size) as a function of the cluster mass in two relevant 
redshift bins: z=0-0.2 (red lines) and z=0.2-0.4 (black 
lines) in a ACDM model. Vertical dashed lines mark 
the [1.9 - 3.8] X lO^^M© mass bin of Fig.(l). Calcula- 
tions are obtained for r]t = 0.3. 

In particular, in agreement with observations and in- 
dipendently from the adopted cosmology, we find that 
20-30% of clusters in the binA can form a RH and that 
only 2-3% of galaxy clusters in the binB host a RH. 
Given the requested values of ry* (Fig.(l)), we find that 
the relatively high occurence of RHs observed in mas- 
sive clusters can be well reproduced by our particle ac- 
celeration model under reasonable conditions, i.e. that 
a fraction of 20-30% of the energy of the turbulent mo- 
tions injected during cluster merger (which corresponds 
to a few percent of the thermal energy) is in the form 
of MS waves. 

III. Radio Halos Statistics and cluster mass 

In the previous Section we have found that the 
probabilty to form a RH has a strong dependence on 
the mass of the host cluster, it goes from few % for 
M < 1.9 X IO^^Mq to 20-30% for M>2 x IO^^Mq, in 
agreement with observations (Fig.(l)). 

Now we calculate the probability to find RHs as a 
function of the cluster mass from our synthetic popu- 
lation of galaxy clusters in a ACDM cosmology: this 
is a crucial point of our model and marks the differ- 
ence with previous studies (e.g., Ensslin & Rottgering 
(2002)). As an example, in Fig. (2) we report this prob- 
ability in two redshift bins z=0-0.2 (black lines) and 
z=0.2-0.4 (red lines). Clusters with mass << IO^^Mq 
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Fig. 3. — Probability to form giant RHs {^^ IMpc 
h'^^ size) with redstiift in tlie mass bin [1.9 — 3.8] x 
10^^ MqH^q in a ACDM cosmology. Calculations are 
performed for j]t =0.3. 

have a negligible probability to form giant RHs; on the 
other hand, such a probability is found to reach ~40% 
for clusters with masses > 3 x IO^^Mq in the 0-0.2 red- 
shift bin. For a comparison with Fig.(l) we also mark 
in Fig.(2) the mass bin [1.9 - 3.8] x IO^^Mq (vertical 
blue dashed Hues). 

Thus, the important finding of these calculations is 
that only massive clusters can host giant RHs {Rh > 
500 kpc h'^Q ) and that the probabilty to form these dif- 
fuse radio sources presents an abrupt increase for clus- 
ters with about M^2 x IO^^Mq (Fig.(l) and Fig. (2)). 
These findings can be simply explained in the frame- 
work of our model. Infact, it can be shown that 
in the C&B model the energy of the turbulence in- 
jected in galaxy clusters is expected to roughly scale 
with the thermal energy of the clusters (Cassano & 
Brunetti 2004; Cassano & Brunetti these proceeding 
Fig(2)). This seems reasonable and immediately im- 
plies that the energy density of the turbulence is an 
increasing function of the mass of the clusters, i.e. 
£t ocT (X M°' (a=0.56-0.67), and thus particle acceler- 
ation is favoured in massive clusters. 

In general, the infall of subclusters through a main 
cluster, which is not very massive, injects turbulence in 
a volume Vt (calculated using Ram Pressure Stripping, 
Fujita, Takizawa, Sarazin 2003; Cassano & Brunetti 
2004) which is found to be smaller than that of giant 
RHs, Vh [Vh = 47ri?|^/3), and thus the efficiency of the 
mechanism is reduced by about a factor of Vt /Vh in the 
case of less massive clusters. On the other hand, major 



mergers between massive subclusters are expected to 
inject turbulence on larger volumes, of the order of Vh, 
and thus the efficiency of the generation of RHs is not 
reduced and this further favour massive objects as the 
parent clusters of RHs. 

More quantitatively, it can be shown (Cassano 
& Brunetti 2004) that the acceleration efhciency % 
(within Vh), triggered by a major merger event scales 
about with X oc (0.75 for M > 3 • W^Mq, 

1.25 for M < IO^^Mq). Since the maximum energy 
of the accelerated electrons is ■^t oc x/{B^ + B'^mb)^ 
where Bcmb = 3.2 ■ (1 -I- z^nG is the strength of the 
equivalent magnetic field of the CMB, and the break 
frequency is Vb oc one has: 

ex Af ^-^-^-^ (1) 

yB +BCMB) 

and consequently massive clusters are statistically 
favourite to have z^b^lO^ MHz (which is the adopted 
condition to define the presence of a RH). 

IV. Evolution of radio halos with redshift 

The probability to form RHs dependes on the combi- 
nation of the energy losses suffered by relativistic elec- 
trons (mainly due to IC losses cx (1 -I- z)^) with the 
acceleration efficiency powered by the turbulence gen- 
erated during cluster mergers (which depends on the 
merger history). 

In this Section we calculate the evolution with red- 
shift of the probability to form RHs in a ACDM cosmol- 
ogy. As an example, in Fig. (3) we report the probabil- 
ity to form a giant RH (~ 700 Kpc in a ACDM cosmol- 
ogy) with redshift in the mass bin: [1.9-3.8] x 10^^ Mq 
for a representative value of rjt (rjt — 0.3). The occur- 
rence of RHs decreases with redshift due to the higher 
IC energy losses. We note however that such a de- 
crease is not dramatic since in a ACDM Universe ma- 
jor mergers develope at slightly higher redshift with 
respect to a EdS Universe. For instance, in the con- 
sidered case the formation rate of RHs is 20-36% at 
relatively low redshift and decreases to 10% at higher 
redhifts (0.3 <z< 0.4). 



V. The Luminosity Functions of Radio Halos 
(RHLFs) 

We have alredy shown that the observed probability 
to find RHs with the cluster mass is well reproduced 
by the C&B model (see Sec. 2). In Cassano & Brunetti 
2004 (see also Cassano & Brunetti these proceedings) 
it has also been shown that the typical synchrotron 
and IC luminosity of RHs can be well reproduced by 
the model assuming that during mergers a few percent 
of the thermal energy of the cluster is in the form of 
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Fig. 4.— RHLFs x P) expected by the C&B 

model. Results are shown for the redshift bins < z < 
0.2 (black curve with points) and 0.2 < z < 0.4 (red 
curve with points). The expected Local RHLF (blue 
curve with points) is also reported together with the 
Local RHLF from Ensslin & Rottgering (2002) (solid 
blue line) for a comparison. 

MS waves {i.e., % > 0.1 — 0.2). Given these promis- 
ing results, in this Section we derive the expected lu- 
minosity functions of giant RHs (RHLFs). First we 
use the probility to form RHs with the cluster's mass 
■^AM (Fig-2) to estimate the mass functions of RHs 
{dNHiz)/dMdV): 



dNniz) ^ dN,i{z) 
dM dV dM dV 



nps X P^li, 



(2) 



where nps = nps{M,z) is the Press & Schechter 
(1974) mass function (we use npg since our model is 
based on Press & Schechter formalism). The RHLF is 
given by: 



dNniz) _ dNniz) j dP^^ 
dV dPiA ~ dM dV / dM ' 



(3) 



In order to derive dPiA/dM in Eq.©, we combine 
the observed correlations between radio power at 1.4 
Ghz [PiAGHz) and bolometric X-ray luminosity {Lx) 
(e.g., Feretti 2003) and between Lx and the virial mass, 
M200 (e.g., Arnaud & Evrard 1999). The used Piaghz- 
M correlation is obtained collecting the data from all 
the known clusters with giant RHs and converting them 
in a ACDM cosmology. In Fig. (4) we report the Lo- 
cal (here calculated for z < 0.05) RHLF (number of 



RHs per Gpc^ as a function of the radio power) ex- 
pected by our model and the expected RHLFs in the 
bins < z < 0.2 and 0.2 < z < 0.4 (hues with points). 
Our RHLFs are compared with the Local {RHLF)eiur 
(blue sohd line) of Ensslin & Rottgering (2002). The 
(RHLF) EhR is obtained combining the X-ray observed 
luminosity function of clusters with the radio luminos- 
ity - X-ray luminosity correlation and assuming that 
a costant fraction frh = 1/3 of galaxy clusters have 
RHs. The most important difference between the two 
luminosity functions is that our RHLF shows a cut- 
off/flattening at low radio powers. We stress that the 
flattening at low powers is a unique feature of particle 
acceleration models since it marks the effect of the de- 
crease of the efficiency of the particles acceleration in 
the case of the less massive galaxy clusters and conse- 
quently the presence of a synchrotron cut-off ui, < 10^ 
MHz. 

Future radio observations (e.g., with LOFAR and 
LWA) should be able to test the presence of such a 
low-power cut-off in the RHLFs and the evolution of 
the RHLFs with redshift. 
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